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Abstract 

An exploratory study for a new determination of the strong coupling constant a s (Mz) 
from existing world data on the proton structure function F% in the Q 2 -range from ~ 5 
to 120 (GeV/c) 2 is presented. The main features of our approach are: 1) the use of 
low-order Nachtmann moments evaluated with a direct contribution from data larger 
than 70% of the total; 2) the inclusion of high-order perturbative effects through the 
soft gluon resummation technique; 3) a direct control over higher-twist effects; and 4) 
the independence from any specific choice of the x-shape of the input parton distribu- 
tions. At next-to-leading order we get a s (M z ) = 0.1209±0.0010 (stat.)± 0.0015 (syst.) 
with a significant dependence upon the order of the moment used. Including soft gluon 
effects we obtain a s (M z ) = 0.1188 ± 0.0010 (stat.) ± 0.0014 (syst.) with a remarkable 
better stability against the order of the moment. Our findings are compared with recent 
DIS determinations of a s (Mz) and directions for future improvements are discussed. 
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1 Introduction 



The coupling constant a s is the only parameter of Quantum ChromoDynamics (QCD) which 
is our present field theory of the strong interactions. The energy (or scale) dependence of 
a s is fixed by the renormalization group equation of QCD, while the value of a s at a given 
scale must be determined from experiments. 

Up to now there are many measurements of a s (/z) at different scales fi, coming mainly 
from analyses of deep inelastic scattering (DIS), e + e~ annihilation and hadron collider 
processes as well as from analyses of the heavy quarkonium system. It is well established 
(cf. Ref. [T]) that the running of a s (/i) is properly predicted by QCD for fi > 1 GeV up to 
4-loop corrections, and therefore it makes sense to extrapolate all the existing measurements 
of a s (n) to a common reference scale, which is conventionally taken to be the rest mass of 
the neutral Z boson, M z = 91.1876 (21) GeV 0. 

The latest world average of a s (Mz), based on results of QCD analyses performed at 
next-to-next-to-leading order (NNLO), is a s {M z ) = 0.1183 ± 0.0027 3J. As far as the 
extraction of a s (Mz) from DIS data is concerned, the current average value at next-to- 
leading order (NLO) is a s (M z ) = 0.119±0.002 (exp.) ±0.003 (th.) where the theoretical 
error includes an estimate of NNLO effects. Recently, two new determinations of a s (Mz) 
have been carried out using the complete NNLO of perturbative QCD, obtaining a s (M z ) = 
0.1166±0.0009 (exp.) ±0.0010 (th.) [4 anda s (M z ) = 0.1143±0.0014 (exp.)±0.0013 {th.) 0. 

The methodology commonly used to extract a s (Mz) form DIS data can be summarized 
into the following steps: 1) fix the perturbative order (and the renormalization scheme) 
of the analysis; 2) choose a parameterization for the input parton distributions at a given 
scale //; 3) fit scaling violations using a s (fJ>) as a free parameter. Finally, one extrapolates 
a s (fi) to a s (Mz). Therefore, the main sources of uncertainties are expected to be: i) the 
effects of high-order perturbative corrections; ii) the dependence upon the specific choice 
of the x-shape of the input parton distributions; and iii) possible effects of higher twists, 
which can mimic scaling violations in a finite Q 2 window. 

The extraction of a s (Mz) from scaling violations of the nucleon structure function F£* 
is characterized by a remarkable abundance of data with reasonably small errors. There is 
however a non-negligible dependence upon the input parton densities, particularly on the 
input used for the gluon distribution. One can speculate that in the non-singlet channel 
the so-called log-slope of has a reduced dependence upon parton distributions. However, 
existing data on the non-singlet part of FJ? , which involves the difference between the proton 
and neutron structure functions, are not very accurate. Alternatively, one can consider the 
structure function F£ p , which is directly a non-singlet quantity; again however the experi- 
mental uncertainties on F% p are typically large (see Ref. [S] for the latest determination of 
a s (M z ) using neutrino DIS data). 

Recently, the strong coupling a s (Mz) has been extracted from the scaling violations 
of truncated moments of the non-singlet part of the structure function F% [7j. The trun- 
cated moments are determined through a parameterization of BCD MS and NMC data 
obtained using a neural network [H], which allows to retain the full experimental informa- 
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tion on errors and correlations. At NLO the authors of Ref. [7j have obtained a s (Mz) = 
0.124 ±[j;88f (exp.) ±g;g^ (t/i.). The method of Refs. 0111 has the nice feature of minimiz- 
ing many theoretical uncertainties; in particular the neural network parameterization of F^ 1 
does not rely on any specific choice for the input parton distributions. In order to avoid 
effects from higher twists, which are contained in the neural network fit, the Q 2 range used 
in Ref. [7] for the extraction of the strong coupling was restricted to Q 2 > 20 (GeV/c) 2 . In 
this way however the sensitivity to scaling violations is reduced and this may generate large 
uncertainties on the extracted value of a s (M z ). 

The aim of this paper is to present an exploratory study for a new extraction of a s (Mz) 
from proton DIS data. The main features of our procedure are: 1) the use of low-order 
Nachtmann moments evaluated with a direct contribution from data larger than 70% of 
the total moment; 2) the inclusion of high-order perturbative effects through the soft 
gluon resummation technique at next-to-leading-log (NLL) accuracy; 3) a direct control 
over higher-twist effects. As a result, our determination of a s (Mz) does not rely on any 
assumed x-shape of the input parton distributions and the Q 2 -range of our analysis can be 
extended from « 5 to ~ 120 (GeV/c) 2 . At next-to-leading order we get a s (M z ) = 0.1209 ± 
0.0010 (stat.) ± 0.0015 (syst.) with a significant dependence upon the order of the moment 
used. Including soft gluon effects we obtain a s (M z ) = 0.1188±0.0010 (stat.) ±0.0014 (syst.) 
with a remarkable better stability against the order n of the moment, at least for n < 8. 

The plan of the paper is as follows. In Section El low-order Nachtmann moments of the 
proton structure function are evaluated in the Q 2 -range 5 < Q 2 (GeV/c) 2 < 120, following 
the procedure of Ref. P . In Section |3] the precise Q 2 -range for a leading twist analysis 
is determined for each moment separately. Then, the pQCD analysis is carried out both 
at NLO and including soft gluon effects at NLL accuracy. In Section IV our findings are 
compared with those of recent DIS determinations of a s [Mz) and directions for future 
improvements are discussed, namely: i) the need for more precise data at least in the 
Q 2 range from 5 to 15 (GeV/c) 2 , which might come from an upgrade of Jefferson Lab to 
12 GeV electron beam; ii) the consideration of more Nachtmann moments in our analysis; 
and iii) the inclusion of all the NNLO effects and of the resummation of large n-logarithms 
at NNLL accuracy. Our conclusions are finally summarized in Section V. 

2 Evaluation of the Nachtmann Moments 

In this Section we review the methodology developed in Ref. [Oj to evaluate the Nachtmann 
moments of the proton structure function starting from world data. There, all the existing 
data in DIS kinematics as well as the new results obtained with the CLAS detector at 
Jefferson Lab in the nucleon resonance production regions, covering a wide Q 2 -range from 
« 0.1 to ~ 5 (GeV/c) 2 , were considered. For sake of completeness of the present paper we 
describe here the methodology of Ref. but we limit ourselves to the world data set at 
Q 2 > 5 {GeV/cf. 

As it is well known ^U|, the Nachtmann moment M n (Q 2 ) of the proton structure function 
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F 2 (x, Q 2 ) is defined as 



rl Cn+1 

M n (Q 2 ) = / dx^F^(x,Q 2 
Jo X 6 



3 + 3(n + l)r + n(n + 2)r 2 



(n + 2)(n + 3) v ' 

where n is the order of the moment, £ = 2x/(l + r) is the Nachtmann variable and 



r = yl + AM 2 x 2 /Q 2 with M being the nucleon mass. We have chosen the Nachtmann 
definition of the moments instead of the Cornwall-Norton one, because with the former all 
the effects arising from target-mass corrections are exactly cancelled out. Hereafter, for sake 
of simplicity, we will omit the superscript "p" and indicate the proton structure function 
simply by F 2 (x, Q 2 ). 

Since moments are fully inclusive quantities, Eq. receives contributions from inelastic 
channels as well as from the elastic peak, namely: 

M n {Q 2 ) = M™{Q 2 ) + M*{Q 2 ) (2) 

with 

M ei (n 2, _ G 2 E (Q 2 ) + r G 2 M (Q 2 ) +1 3 + 3(n + l)r et + n(n + 2)r 2 el 

M " W) " \T~r * d (n + 2)(n + 3) {6) 

where Ge(Q 2 {Gm{Q 2 )) is the charge (magnetic) Sachs proton form factor, £ el = 2/(1 +r e i), 
r el = Jl + 1/r and r = Q 2 /AM 2 . 



The evaluation of the inelastic moment M™(Q 2 ) involves the computation at fixed Q 2 
of an integral over the whole range of values of x up to the inelastic pion threshold. For 
this purpose world data on the structure function F 2 from Refs. [TTJ H21 EH UH EEH1 HE1 HZj 
and data on the inelastic cross section ^HJ EJ] are considered in order to reach an adequate 
kinematical coverage. The integral over x is performed numerically, using the standard 
trapezoidal method TRAPER [20] ■ Data from Ref . |2I] are not included in the analysis due 
to their inconsistency with other data sets as explained in detail in Ref. [22] > and data from 
Ref. j2Sl 123! are n °t considered due to the large experimental uncertainties. 

First of all the Q 2 -range from 5 to 120 (GeV/c) 2 is divided into bins having a width equal 
to AQ 2 /Q 2 = 5%. Bins do not follow one another, but they are located only where data 
exist. Nevertheless, different bins do not overlap. Then world data are sorted within each 
Q 2 bin and shifted to the central bin value Qq, using the fit from Ref. [23], which consists of 
the SMC parameterization of Ref. [21 at W > 2.5 GeV and the SLAC fit from Ref. [IE] 
for W < 2.5 GeV, namely 



F 2 {x,Ql)= Ft(x,Q 2 ). (4) 



F 2 (x,Q 
F{*(x,Q 

Then the difference between real and bin-centered data 



5^\x,Ql)=F 2 {x,Q 2 ) 



(5) 
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is added (in quadrature) to the systematic errors of F 2 . 

In order to have a data set dense in x, which reduces the error in the numerical integration, 
we perform an interpolation, at each fixed Qq, when two contiguous experimental data points 
differ by more than V. The value of V is not unique and depends on kinematics. We divide 
the whole range of x (or W) in three intervals: resonance region (W < 1.8 GeV), normal 
DIS region (W > 1.8 GeV and x < 0.06) and very \ow-x region (x < 0.06). In the 
resonance regions, where the structure function exhibits strong variations, we put V equal 
to 0.01. Above resonances where F 2 is smooth, V = 0.1, and finally in the low x region we 
set V = 0.005. Changing these V values by as much as a factor of two produces changes in 
the moments that are much smaller than the final systematic errors reported in Table H 

To fill the gap between two contiguous data points, x a and Xb, which are distant each other 
more than V, we use the interpolating function F 2 mt {x 1 Q 2 ) defined as the parameterization 
from Ref. [22] normalized to the data on both edges of the interpolating range. Assuming 
that the x-shape of F 2 nt (x,Q 2 ) is correct, one has 

F™ t (x,Q 2 ) = Il(Q 2 )F 2 fu (x,Q 2 ), (6) 

where the normalization factor II (Q 2 ) is defined as the weighted average evaluated using all 
experimental points located within the interval A around x a or xj,, namely 



n(Q 2 ) = 4(Q 2 ) E r .... (7) 



5p 2 at (x h Q 2 

where 5 s p^{xj,Q 2 ) is the statistical error relative to F 2 lt and 

\ -1/2 

\xi-x ab \<A 



Sn(q 2 ) = E 7 T2 

[Sg at (x t ,Q 2 ) 

is the statistical uncertainty of the normalization. Therefore, the statistical error of the 
moments calculated according the trapezoidal rule [20] was increased by adding the linearly 
correlated contribution from each interpolation interval as follows: 

S n (Qo) - 6n(Q ) J dx— F 2 (x, Q ) {n + 2){n + 3) • ( 9 ) 

Since we average the ratio F 2 (xi,Q 2 )/F2 lt (xi,Q 2 ), the value of A is not affected by 
the resonance structures in F 2 , and we fix it by the requirement to have more than two 
experimental points in most cases, obtaining A = 0.01 in the resonance and in the very 
low-x regions, and A = 0.1 in the DIS region. In Fig. Q]we show how this interpolation is 
applied: the thin continuous line represents the original function F 2 lt (x,Q 2 ) and the heavy 
continuous line represents the result of the interpolation F 2 nt (x,Q 2 ). We have also checked 
that the moments do not show any significant dependence on the specific value of A. 
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Figure 1: Example of the interpolating procedure. The meaning of both the curves and the symbols 
is described in the text. 

To fill the gap between the last experimental point and one of the integration limits 
(x a = or Xb = 1) we perform an extrapolation at each fixed Q 2 using F 2 ^*(x, Q 2 ) including 
its systematic error given in Refs. (25J CHI I2H1 122] ■ 

The total systematic error consists of genuine uncertainties in the data, as given in 
Refs. [TTJ H21 E3 EH US1 EE 113 UH EH], and uncertainties in the evaluation procedure. To 
estimate the first type of error one has to take into account that many data sets have been 
measured in different laboratories and with different detectors. We assume that different 
experiments are independent and therefore only systematic errors within one data set are 
correlated. Thus, an upper limit for the contribution of the systematic error from each data 
set is evaluated in the following way: 

• we first apply a simultaneous shift to all experimental points in a data set by an amount 
equal to their systematic error; 

• then the inelastic n-th moments obtained using these distorted data, M™ od (Q 2 ), are 
compared with the original ones M™(Q 2 ) evaluated with no systematic shifts; 

• finally the deviations for each data set are summed in quadrature as independent 
values. 

Explicitly, one has 

~N~ S 

E [Mr d (Q 2 ) - m™(Q 2 )] 2 (io) 




M™{Q 2 ) \ 
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where Ns is the number of available data sets. Then the resulting error is summed in 
quadrature to 5™ rm (Q 2 ). 

The second type of error is related to bin centering, interpolation and extrapolation. The 
systematic uncertainty due to bin centering is estimated as follows: 

5 c n {Q 2 ) =Y,K n (x l ,Q 2 ) Wl {Q 2 ) 6 c ^\x t ,Q 2 ) (11) 

i 

where according the Nachtmann moment definition and the trapezoidal integration rule one 
has 

K n ( Xi ,Q 2 ) = g+ 1 3 + 3(n + l)r i + n(n + 2)rf 



x\ (n + 2)(n + 3) 

Wi(Q 2 ) = (x i+1 - ari_i)/2 . (12) 

The relative systematic error of the interpolation is estimated as the possible change of 
the fitting function slope in the interpolation interval and it is evaluated as a difference in 
the normalizations at different edges: 



1 |a ^ <A F 2 ( Xl ,Q 2 ) 1 '^ <A F 2 (x v Q 2 



(13) 



where Ni and Nj are the number of points used to evaluate the sums. Since the relative 
function F 2 j 1 F 2 lt is a very smooth function of x in DIS kinematics, the linear approximation 
can provide a good estimate on a limited x-interval (smaller than V). Thus, the error given 
by Eq. ()13|) accounts for such a linear mismatch between the fitting function and the data 
on the interpolation interval. Therefore, the systematic error introduced in the moments by 
the interpolation procedure can be estimated as 



b 



&i(ri*\ x (ri*\ 1a jrfitf n*\ 3 + 3(n+l)r + n(n + 2)r 2 
S n (Q ) = 5 S (Q ) J dx— F( (x, Q ) (n + 2)(n + 3) ■ (14) 

The systematic errors obtained by all these procedures are summed in quadrature: 



$n(Q 2 ) = vtTO 2 )] 2 + mQ 2 )} 2 + mQ 2 )? ■ (is) 

In order to study the systematic error due to extrapolation in unmeasured regions we have 
performed a test on a possible dependence upon the functional form of the fitting function 
by comparing the moments presented here with those obtained using the neural network 
parameterization of Ref. 0. The difference is expected to be significant only for the second 
moment M 2 (Q 2 ), due to the presence of a factor x n ~ 2 in the integrand of the Nachtmann 
moments. Any way we have checked the difference for all the moments considered in this 
work. Such a difference turns out to be smaller than 5^{Q 2 ) given by Eq. (JT5j) [see Fig. |2] 
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Figure 2: Errors of the inelastic Nachtmann moment M^ijQ 2 ) in percentage: empty circles repre- 
sent statistical errors; empty crosses show the systematic error obtained in Eq. il5\) : the difference 
between inelastic moments extracted using the two different F2 parameterizations of Ref. \23$ and 
Ref. is shown by stars. 



both for n = 2 and n > 4. The difference is added to 5^{Q 2 ) in quadrature to evaluate the 
total systematic error of the n-th moment. 

According to Eq. (J2J the contribution from the proton elastic peak should be added to 
the inelastic moments so far obtained. The Q 2 -behavior of the proton elastic form factor 
is parameterize as in Ref. |2"T| . modified accordingly to the recent data on Ge/Gm [28J . 
The uncertainty on the form factors is taken to be equal to 3% according to the analysis of 
Ref. [29], and it is added quadratically to both the statistic and the systematic errors. The 
elastic contribution M^{Q 2 ) turns out to be a quite small correction for Q 2 > 5 (GeV/c) 2 . 
Our final results for the total (inelastic + elastic) moments with n = 2, 4, 6, 8 including both 
statistic and systematic errors, are shown in Fig. El and reported in Table d It can cleary 
be seen that scaling violations are quite small in the second moment M 2 (Q 2 ), while they 
increase significantly as the order n increases. Consequently, any extraction of a s (Mz) from 
the second moment is expected to be plagued by very large uncertainties, as we have explicitly 
verified. It is for this reason that in the next Section the QCD analysis of scaling violations 
will be restricted to n > 4. Note also that the amount of the direct experimental contribution 
to M n (Q 2 ) reduces between 50% and 70% at large Q 2 \ in particular: at Q 2 > 75 (GeV/c) 2 
for n = 4 and n — 6, but already at Q 2 ~ 20 (GeV/c) 2 for n — 8. Finally the systematic 
uncertainties increase significantly as Q 2 increases. 

Before closing this Section, we want to stress two important points: 
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0.20 I 




Q 2 (GeV/c) 2 

Figure 3: Total (inelastic + elastic) Nachtmann moments M n (Q 2 ) [Eq. extracted from the 
proton world data in the range 5 < Q 2 < 120 (GeV/c) 2 . Open dots, squares, diamonds and triangles 
correspond to n = 2, 4, 6 and 8, respectively. Both statistical and systematic errors are reported; 
the former are those slightly shifted to the right side of the markers. 

• from Fig. El it can clearly be seen that for n > 4 the Q 2 -slope of the moments changes 
significantly around Q 2 ~ 20 -f- 30 (GeV/c) 2 . In our opinion this should be related to 
the crossing of the 6-quark threshold, Q 2 , which here after we fix at Q 2 = 25 (GeV/c) 2 . 
The Q 2 -region between » 5 and » 20 (GeV/c) 2 turns out to be the most sensitive one 
to scaling violations and, therefore, to the effects of high-order radiative corrections, as 
already shown in Ref. (SO] (see in particular Fig. 1 of Ref. [30 ). Therefore we consider 
crucial to include such a region for a precise determination of a s (M z ), provided possible 
power correction effects are taken care of; 

• thanks to the direct use of data as well as of parameterizations like the one of Ref. [B] 
our evaluation of the proton Nachtmann moments is independent of any specific choice 
of the x-shape of input parton distributions. 



3 QCD analysis of proton moments with n > 4 

The Q 2 -behavior of the moments M n (Q 2 ) can be analyzed adopting the powerful tool of the 
Operator Product Expansion (OPE). Our extracted Nachtmann moments can be expanded 
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Table 1: Total (inelastic + elastic) Nachtmann moments M n (Q 2 ) for n = 2,4, 6 and 8 evaluated 
in the interval 5 < Q 2 < 120 (GeV/c) 2 . The results, reported together with their statistical and 
systematic errors, are labeled with an asterisk when the contribution to the integral by the experi- 
mental data is between 50% and 70% of the total; all other results are obtained with more than 70% 
direct contribution from data. Adapted from Ref. [9]. 



Q l(GeV/cy\ 


AfttQ^xlO- 1 




M 6 (Q^)xlO- J 




5.967 




1.860 ± 0.015 ± 0.116 


4.995 ± 0.046 ± 0.553 




7.268 




1.768 ± 0.012 ± 0.044 






7.646 






4.552 ± 0.045 ± 0.098 


1.714 ± 0.021 ± 0.040 


8.027 






4.429 ± 0.027 ± 0.135 


1.647 ± 0.011 ± 0.064 


8.434 




1.668 ± 0.014 ± 0.084 


4.350 ± 0.033 ± 0.109 


1.633 ± 0.013 ± 0.042 


8.857 


1.724 ± 0.015 ± 0.041 


1.658 ± 0.019 ± 0.027 


4.215 ± 0.042 ± 0.109 


1.562 ± 0.014 ± 0.057 


9.781 




1.662 ± 0.010 ± 0.061 


4.205 ± 0.034 ± 0.146 


1.546 ± 0.017 ± 0.067 


10.267 


1.753 ± 0.015 ± 0.052 


1.629 ± 0.019 ± 0.031 


4.079 ± 0.035 ± 0.095 


1.486 ± 0.011 ± 0.037 


10.793 






4.039 ± 0.106 ± 0.761 


1.471 ± 0.018 ± 0.110 


11.345 


1.732 ± 0.016 ± 0.041 


1.578 ± 0.018 ± 0.035 


3.896 ± 0.041 ± 0.140 


1.402 ± 0.018 ± 0.079 


11.939 


1.759 ± 0.008 ± 0.042 


1.600 ± 0.013 ± 0.031 


3.946 ± 0.040 ± 0.111 


1.418 ± 0.018 ± 0.057 


13.185 




1.528 ± 0.016 ± 0.021 


3.705 ± 0.029 ± 0.067 


1.323 ± 0.009 ± 0.034 


15.310 


1.686 ± 0.014 ± 0.074 


1.473 ± 0.019 ± 0.032 


3.547 ± 0.044 ± 0.133 




16.902 


1.718 ± 0.010 ± 0.051 


1.451 ± 0.017 ± 0.025 


3.401 ± 0.058 ± 0.073 


1.183 ± 0.025 ± 0.028 


18.697 


1.679 ± 0.033 ± 0.097 


1.402 ± 0.013 ± 0.027 


3.281 ± 0.039 ± 0.088 


1.144 ± 0.015 ± 0.035 


19.629 








♦1.111 ± 0.022 ± 0.047 


21.625 


1.677 ± 0.008 ± 0.041 








24.192 


n.711 ± 0.007 ± 0.114 


1.385 ± 0.008 ± 0.024 




♦1.077 ± 0.014 ± 0.043 


26.599 


♦1.665 ± 0.062 ± 0.135 








28.192 


♦1.702 ± 0.009 ± 0.140 


1.344 ± 0.007 ± 0.037 




♦1.012 ± 0.010 ± 0.081 


31.858 


♦1.703 ± 0.010 ± 0.096 








36.750 


♦1.696 ± 0.013 ± 0.111 


1.314 ± 0.009 ± 0.057 


2.971 ± 0.027 ± 0.313 


♦1.003 ± 0.014 ± 0.191 


44.000 


♦1.681 ± 0.013 ± 0.085 








49.750 


♦1.658 ± 0.019 ± 0.101 








57.000 


*1.694 ± 0.017 ± 0.170 








64.884 


♦1.636 ± 0.043 ± 0.114 


1.222 ± 0.053 ± 0.044 


2.708 ± 0.082 ± 0.193 


♦0.895 ± 0.016 ± 0.116 


75.000 




*1.206 ± 0.008 ± 0.025 


♦2.651 ± 0.024 ± 0.150 




88.000 


♦1.669 ± 0.088 ± 0.075 


♦1.199 ± 0.038 ± 0.035 


♦2.630 ± 0.057 ± 0.202 




99.000 




♦1.179 ± 0.012 ± 0.034 


♦2.568 ± 0.029 ± 0.228 




117.75 


♦1.566 ± 0.051 ± 0.088 


♦1.128 ± 0.050 ± 0.052 


♦2.457 ± 0.162 ± 0.261 





as a series of twist terms 

M n (Q 2 ) = r ]n (Q 2 ) + RT n (Q 2 ) (16) 

where i] n (Q 2 ) is the leading twist moment of order n and HT n (Q 2 ) stands for the total 
higher-twist contribution, which is power suppressed. 

The twist- 2 moments f] n {Q 2 ) are generally given by the sum of a singlet and non-singlet 
(NS) terms. The corresponding anomalous dimensions are fixed by pQCD evolution and 
therefore one is left with a total of three unknown parameters, namely the values of the 
gluon, singlet and NS quark moments at a given scale /i 2 . However, the decoupling in the 
pQCD evolution of the singlet quark and gluon densities at large x allows to consider a 
pure NS evolution for n > 4 (cf. [25 J. This means that, besides the quantity a s (fi 2 ), we 
have only one twist-2 parameter for n > 4, namely the value of the twist- 2 moment rj n {ix 2 ) 
at the chosen scale jj 2 . The NS approximation is expected to be quite reasonable for an 
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exploratory study, like the present one. The inclusion of the full singlet and NS evolutions 
for the moments with n > 4 is postponed to a future work. 

As far as the order of perturbation theory is concerned, the NLO approximation is not 
reliable for moments with order n > 4 [30!, which are dominated by kinematical regions 
corresponding to large values of x. Therefore we estimate radiative corrections beyond the 
NLO according to the soft-gluon resummation (SGR) technique developed in Refs. pTr] l32]. 
We stress that SGR, unlike renormalon models, does not introduce any further parameter 
in the description of the leading twist. With respect to Ref. [20], where SGR was considered 
for the quark coefficient function only, we consistently add in this work the resummation of 
large n- logarithms appearing also in the one-loop and two-loop NS anomalous dimensions' 3 . 
Explicitly, for n > 4 one has 

3 ia s (Q 2 " 



Vn(Q' 



+ 



A n [a s (Q 2 ) 
a s (Q 2 



4tt 



-R 



NS 



1 + 



47T 



2C 



{NLO) 
DIS 



+ A 7 



(1,NS)\ 
DIS ) 



G n (Q 2 ) 



(17) 



where A n is a constant absorbing in its definition the n-th moment of the leading twist at 
the scale fi 2 , and is the one-loop NS anomalous dimension. In Eq. (fT7|) the quantity 
R^ s is given by 



R 



NS 



n (NLO) n {NLO) n (NLO) 
°n ^DIS ^n,LOG 



+ A 7; 



(l,NS) 



A7 



(1,NS) 
DIS 



- A 7 , 



(1,NS) 
LOG 



(18) 



where 



A7, 



(1,NS) 



(1,NS) _ Pi NS 
in a 'n 



(19) 

11 - 2A7/3, & = 102 - 



with , y n 1 ' NS ^ being the two-loop NS anomalous dimension, /3 = 
38Nf/3 and Nf the number of active quark flavors at the scale Q 2 . 

In Eq. ()18|) C n NL °^ is the NLO part of the quark coefficient function, that in the MS 
scheme (which is adopted throughout this paper) reads as 



C^ = C F \S 1 (n) 



Si(n) + ^ 1 

2 nln + 1, 



S 2 (n) + f + + (20) 
2n n + 1 n z 2 



where Cp = (N 2 — 1) / (2N C ) and Sk{n) = S"=i l/j fc - F° r large n (corresponding to the large- 



x 



region) the coefficient C n NL °^ is logarithmically divergent; indeed, since Si(n) = 7B+log(n) + 
0(l/n), where 7 s = 0.577216 is the Euler-Mascheroni constant, and Sz(n) = iT 2 /6 + 0(l/n), 
one gets 



c (NLO) 



C^ + C^ + Otl/n) 



i(NLO) 



C 



(NLO) 
DIS 



c f 



7T 



n {NLO) 
^n,LOG 



C F ln(n) 



ln(n) + 2 lE + - 



(21) 



3 Such an addition turns out to have a quite small effect (+0.0003) on our final determination of a s {Mz 
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In case of the quantity A74 • ^ one obtains 

(1,NS) a 



A7 



n 




A„ 


(1,N S) 






DIS 




A) 




+ 


C A 




+ 


N f 


(1,NS) 




c F 


n,LOG 




A> 



(l,NS) 



c h 



2vr 2 + 32S(oo) -4S 3 (oo) 



22 ~ 17 

f - 2 " 165(«») - £ 



4tt 2 



+ 



—I 



ln(n) , 



(22) 



where C A = iV c , 5(oo) = E^i(-)^i(j)/j 2 = -0.751286, S 3 (oc) = 1.202057 and K = 
C A (67/18 - tt 2 /6) - 5^/9. 

In Eq. ()17|) the function G n (Q 2 ) is the key quantity of the soft gluon resummation. At 
NLL accuracy it reads as 



G n (Q 2 ) = ln(n) Gi(X n ) + G 2 (X n ) + O a k s \n k ~\n) 



(23) 



where A„ = fi a s (Q 2 ) ln(n)/4?r and 



Gi(A) 



G 2 (A) 



C f 



A)A 



[A + (1 - A)ln(l - A)] 



47 g + 3 
-C F — ln(l - A) 

Po 

+^ln(l 

Po 



A) 



CvSn(l 
Po 



-ln(l - A) 



A) 



(24) 



with K = N c (67/18 - tt 2 /6) - 5N f /9. Note that the function G 2 (A) is divergent for 
A — *■ 1; this means that at large n (i.e. large x) the soft gluon resummation cannot be 
extended to arbitrarily low values of Q 2 . Moreover, the sub-leading terms 0[agln fc_1 (n)] 
become more and more important as A n increases. Therefore, for a safe use of the SGR 
technique at NLL accuracy it is essential to check that X n is small enough. Later on, we 
will impose the condition that the Q 2 -ranges of our analysis should correspond to the cut 
X n < 0.3; this appears to be a safe upper limit at least to avoid significant effects from 
NNLL terms in Eq. (J23j) . as we have checked using the NNLL order expression for the 
function G n (Q 2 ) reported in Ref. [3*3] . It is straightforward to see that in the limit A„ << 1 

one has G n (Q 2 ) -»■ a s (Q 2 ) [2 C^og + ^Infoa]/^, so that Ec i- d reduces to the well- 
known NLO approximation. This implies that adopting the two-loop approximation for the 
running coupling constant a s (Q 2 ) the twist- 2 expression (fTTj) contains all the NLO effects 
and the resummation of large n- logarithms beyond the NLO. 
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To sum up, the leading twist (JT7j) contains two unknown parameters: the coefficient A n 
and the value of the coupling constant a s (fi 2 ) at the scale [i 2 . Alternatively, as another 
parameterization for a s (Q 2 ), the so-called QCD scale parameter A is introduced in a such 
a way that at two- loops the running of the coupling constant a s (Q 2 ) can be written as 

a (Q 2 ) - ^ jl - A lnKQVA 2 )] ! (25) 

as{Q } ~ (3 HQ 2 /A 2 ) V ft HQ 2 /A 2 ) ) ' {Zb) 



Actually the parameter A depends on the renormalization scheme (the MS one in this paper) 
and on the number of active quark flavors. Therefore, more properly, in Eq. (|25j) A should 
be labelled as A^- . The constraints of continuity of a s {Q 2 ) at the various quark thresholds 

impose relations between the parameters A^^ for the various values of Nf, leaving only one 
of them as a free parameter to be determined from the data. 

As far as power corrections are concerned, several higher-twist (HT) operators exist and 
mix under the renormalization-group equations. Such mixings are rather involved and the 
number of mixing operators increases rapidly with the order n of the moment. A complete 
calculation of the HT anomalous dimensions for n > 2 is not yet available, and therefore a 
phenomenological ansatz is usually adopted for BT n (Q 2 ) in the literature (cf. Refs. [213 HUB 
|5]), namely 



KT n (Q 2 ) = ^ [a s (Q 2 )\ + ^ [a s (Q 2 )\ + ... , (26) 

where the logarithmic pQCD evolution of the twist-r contribution is accounted for by the 
term [a s (Q 2 )\ 1n with an effective anomalous dimension 7^ and the parameter rep- 
resents the overall strength of the twist-r term. Finally, in Eq. (|26| the ellipses stand for 
contributions of twists higher than r = 6. 

The twist analysis of Ref. [§] includes the new precise data set obtained with the detector 
CLAS at Jefferson Lab, which dominates the Q 2 -region below w 5 (GeV/c) 2 . In Ref. jH], 
where all the HT parameters appearing in Eq. ()26j) were determined from data fitting, it 
was found that power corrections are relevant in the few (GeV/c) 2 region, and for Q 2 > 
1 (GeV/c) 2 they can be described fairly well by the first two terms of the HT expansion 
(J2HJ). The Q 2 -behavior of the ratio between the total HT term ()26|) and the total moments 
M n (Q 2 ), as determined in Ref. 0, is reported in Fig. E|for Q 2 > 5 (GeV/c) 2 . 

In Ref. jH] the value of a s (M z ) was not determined from the data, but it was kept fixed 
close to the world average, namely a s [Mz) = 0.1180. The reason for that is the occurrence 
of an (anti) correlation between the value of a s (Mz) and the total HT strength. Such a 
correlation does not allow to extract a s (Mz) and the HT simultaneously from our proton 
data with good precision. This point was already addressed in Ref. [30j, where the impact 
of the specific value of a s [Mz) on the extracted higher twists was investigated. There it 
was shown that the strength of the total higher-twist decreases as a s (Mz) increases without 
changing the statistical significance of the fit; indeed, the calculated x 2 turned out to be 
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tf (GeV/c) 2 

Figure 4: i?a£zo in percentage between the total HT contribution HT n (Q 2 ) [Eq. \2b)) l and the 
moments M n (Q 2 ) for Q 2 > 5 (GeV/c) 2 . The HT contribution has been obtained in Ref. -9f using 
proton world data in the Q 2 -range 1 < Q 2 (GeV/c) 2 < 100. The value of a s (Mz) is kept fixed at 
a s (Mz) = 0.1180. Open squares, diamonds and triangles correspond to n = 4, 6 and 8, respectively. 
Error bars represent the uncertainties resulting from the fitting procedure of Ref. Adapted from 

RefiW- 

almost insensitive to the specific combination of the values of a s (Mz) and of the total HT 
strength. 

The previous considerations clearly suggest that for the extraction of a s (Mz) from the 
moments M n {Q 2 ) a pure leading- twist analysis is necessary, but the lower limit Q^in °f the 
range of the analysis should be chosen in such a way that perturbative effects are maximized 
for each order n separately, while HT effects can be still considered under control. Both 
goals can be achieved as follows: 

• since we work within the NLL approximation for the resummation of large n-logarithms, 
we have to impose the condition that NNLL effects are negligible, which means that 
we should keep the parameter A n = j3 a s (Q 2 ) m(n)/47T, appearing in Eq. (j23J), small 
enough. Effects of the NNLL order on the SGR quantity G n {Q 2 ) have been calcu- 
lated in Ref. [S3] and we have checked that the cut A n < 0.3 makes those NNLL 
effects negligible for n = 4, 6 and 8. We point out however that our present estimate 
of NNLL effects is not complete. Effects of subleading logs have been recently investi- 
gated in Ref. [HI], where an upper bound A n ~ 0.4 was assumed at NNLL accuracy. At 
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a s (Mz) = 0.1180 our upper bound A n ~ 0.3 corresponds to Q 2 ~ 5, 7 and 11 (GeV/c) 2 
for n = 4, 6 and 8, respectively; 

• from Table ^ the average value of the systematic errors on our moments in the Q 2 - 
range from » 5 to ~ 20 (GeV/c) 2 turns out to be f=a 3% at n = 4, 6 and ~ 4% at 
n = 8. Looking now at Fig. H]it can be seen that such percentage values are reached 
by the total HT contribution when Q 2 9, 10 and 16 (GeV/c) 2 for n = 4, 6 and 8, 
respectively. 

Therefore, for the lower limit Q 2 min we will consider the following ranges of values in units 
of [(GeV/c) 2 }: 5 < Q 2 mm < 9 for n = 4; 7 < Q 2 mm < 10 for n = 6; and 11 < Q 2 mm < 16 
for n = 8. At the same time, in order to be as much conservative as possible, the total HT 
strength shown in Fig. 0] is added (in quadrature) to the systematic error of our moments. 
Finally, the upper limit of the Q 2 -range of our analysis is taken to be equal to 120 (GeV/c) 2 
for each n according to the data set of Table ^ 

Using a ^-minimization procedure, we have fitted the moments M n (Q 2 ) neglecting the 
HT contribution HT n {Q 2 ) in Eq. (JTfij) and adopting the SGR formula (fTTj) for the leading- 
twist term rj n (Q 2 ). The results obtained at n = 4, 6 and 8 for various values of Q 2 min are 
reported in Table El There, the columns labelled "stat." represent the uncertainties due 
to the fitting procedure and driven by the statistical errors of our moments. The columns 
labelled "syst." represent our estimate of the systematic uncertainties on the extracted values 
of a s (Mz), obtained simply by adding (subtracting) to the central values of the moments 
their own systematical errors and by repeating the ^-minimization procedure. In this way 
an upper (lower) value of a 3 (Mz) is generated by the systematic errors of the moments. 

From Table |21 it can be seen that for n = 4 and n = 6 all the determinations of a s (Mz) 
are consistent with each other within the statistical error. The systematic uncertainties are 
quite limited and not larger than the statistical ones. The same does not hold in case of 
n — 8. Indeed, using Q 2 min = 11,12 and 13 (GeV/c) 2 we still extract values of a s (Mz) 
almost consistent with each other, but larger than those obtained at n = 4, 6. Moreover, for 
Qmin = ^ (GeV/c) 2 the determination of a s (M z ) drops down to a very low value. Finally, 
the systematic uncertainties become quite large at n = 8. Note that: i) the number of data 
points is the smallest one at n = 8; ii) the amount of the direct experimental contribution 
to M$(Q 2 ) reduces between 50% and 70% already at Q 2 ~ 20 (GeV/c) 2 with large systematic 
uncertainties (see Table HJ). Thus, the results obtained at n = 8 may be explained, at least 
partially, by the limited quality of the data, but they might also represent the signal for the 
need of including higher-order perturbative effects. 

For each moment we have separately calculated the (weighted) average of the extracted 
values of a s {Mz) over the various values of Q 2 min - The results are collected in Table El The 
central values for n = 4, 6 are in good agreement with each other and, surprisingly, also the 
one at n = 8. However, we do not consider significant such a coincidence because of the 
spread of the results reported in Tabled at n — 8. Then, we have calculated the (weighted) 
average of the results obtained for the various values of n, and we have reported them again 
in Table El 
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Table 2: Values of a s (Mz) extracted from a pure leading-twist analysis of the moments M n (Q 2 ), 
using the SGR formula for various values of the lower limit Q min of the analysis range, 

as described in the text. The values of Q min are in (GeV/c) 2 , while the upper limit is fixed at 
Qmax = 120 (GeV/c) 2 . The columns labelled "stat." and "syst." represent the uncertainties due to 
the fitting procedure and to the systematic errors of the moments M n (Q 2 ), respectively (see text). 



n = 4 



O' 2 ■ 


N° of data 


a s {M z ) 


stat. 


syst. 


5 


20 


0.1198 


0.0027 


0.0009 


7 


19 


0.1191 


0.0030 


0.0005 


8 


18 


0.1179 


0.0035 


0.0009 


9 


16 


0.1183 


0.0037 


0.0006 



n = 6 



O' 2 ■ 


N° of data 


a s (M z ) 


stat. 


syst. 


7 


19 


0.1197 


0.0022 


0.0011 


8 


18 


0.1194 


0.0022 


0.0013 


9 


15 


0.1181 


0.0027 


0.0025 


10 


14 


0.1169 


0.0030 


0.0032 



n = 8 



O 2 ■ 


N° of data 


a s (M z ) 


stat. 


syst. 


11 


10 


0.1214 


0.0030 


0.0050 


12 


9 


0.1213 


0.0027 


0.0078 


13 


8 


0.1196 


0.0034 


0.0088 


16 


7 


0.1085 


0.0039 


0.0105 



As our final estimate of a s (M z ) we exclude the determination at n = 8 and consider only 
the (weighted) average of the results at n = 4 and n — 6, leading to 

a s (M z ) = 0.1188 ± 0.0010 (stat.) ± 0.0014 (syst.) (27) 

We have checked that a change in the 6-quark threshold Q 2 between « 20 and ~ 30 (GeV/c) 2 
introduces an uncertainty of ~ ±0.0005, which can be neglected in comparison with the sys- 
tematic error reported in Eq. (|27p. It should be mentioned that the uncertainty in the extrac- 
tion of a s (M z ) may be affected also by the renormalization scale dependence of Eq. (fTTj) . 
Such a dependence is known in literature j33 and it is a common practice to investigate 
its impact on the extracted value of a s (M z ), obtaining in this way an estimate of higher- 
order perturbative effects. However, due to the exploratory nature of our present work, we 
postpone such a (mandatory) investigation to a future work. 

In order to disentangle the effects of the soft-gluon resummation we have repeated our 
calculations adopting the NLO approximation for the leading twist i] n (Q 2 ). The results 
are reported in the Tables It can be seen that our NLO determination of a s (M z ) 
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Table 3: Weighted averages of the values of a s (Mz) reported in Table\^for each moment M n (Q 2 ) 
separately. The last two rows represent the final (weighted) averages of a s (M z ) taken over the 
various moments considered. 



moment 


a s {M z ) 


stat. 


syst. 


M 4 


0.1189 


0.0016 


0.0007 


M 6 


0.1188 


0.0012 


0.0018 


M 8 


0.1189 


0.0016 


0.0077 






average 


0.1189 


0.0008 


0.0031 


Mi + M 6 


0.1188 


0.0010 


0.0014 



exhibit features similar to the ones observed for the SGR results of Tables I2H2J including 
the disturbing strong dependence on Q 2 min at n = 8, but with an important difference: the 
determinations of a s {Mz) for the various values of n are not consistent with each other and 
a remarkable dependence on the order of the moment is clearly visible in Table E] and in 
Fig. El where the results obtained at NLO and within the SGR are compared directly. 

In Fig. El our final determination of a s (Mz) [see Eq. (12 7j) ] is compared with the total 
and DIS world averages, taken from Ref . [3] . It can be seen that the comparison is positive 
and our central value agrees very well with the world averages within the quoted uncertain- 
ties. In the same figure we have included also the results of the recent NNLO analyses of 
Refs. Jl] and [5j. The latter appears to be significantly lower than our determination (|27|) . 
while the result of Ref. |3J is still consistent with our result but at the limits of the quoted 
uncertainties. We expect that an inclusion of all the NNLO effects with a consistent resum- 
mation of large n-logarithms at NNLL accuracy will bring our result closer to the one of 
Ref. [lj and possibly also to the findings of Ref. j^j. We point out however that the result 
of Ref. [5] has been obtained with a procedure which depends on the specific choice of the 
input parton distributions, while our procedure and the one of Ref. jl] are totally free from 
such a dependence. 

The results presented in our exploratory study show that the new approach appears to 
work quite well and it may represent an interesting tool for providing a precise determina- 
tion of a s (Mz) from DIS data. Clearly our approach should be improved both from the 
experimental and the theoretical sides. 

As far as data are concerned, we need a more precise determination of the Nachtmann 
moments in the Q 2 range from w 5 to ~ 20 (GeV/c) 2 , where the largest sensitivity to scaling 
violations is expected. In this respect we point out that thanks to its large acceptance the 
CLAS detector at Jefferson Lab has already demonstrated its capability to determine 
the cross section in a wide two-dimensional range of values of x and Q 2 , allowing a direct 
integration of the data at fixed Q 2 over the whole significant x-range for the determination of 
the Nachtmann moments with n > 2. Presently with the 6 GeV electron beam, operating at 
Jefferson Lab, the upper limit in Q 2 is around 5 (GeV/c) 2 . The proposed upgrade to 12 GeV 
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Table 4: The same as in Tabled but adopting the NLO approximation for the leading twist. 



n = 4 



O 2 ■ 


N° of data 


a 8 {M z ) 


stat. 


syst. 


5 


20 


0.1201 


0.0027 


0.0008 


7 


19 


0.1194 


0.0029 


0.0004 


8 


18 


0.1182 


0.0034 


0.0004 


9 


16 


0.1186 


0.0037 


0.0006 



n = 6 



O 2 ■ 


N° of data 


a 8 {M z ) 


stat. 


syst. 


7 


19 


0.1232 


0.0024 


0.0013 


8 


18 


0.1228 


0.0025 


0.0014 


9 


15 


0.1214 


0.0029 


0.0028 


10 


14 


0.1200 


0.0032 


0.0035 



n = 8 



O' 2 ■ 


N° of data 


a s (M z ) 


stat. 


syst. 


11 


10 


0.1276 


0.0030 


0.0060 


12 


9 


0.1273 


0.0034 


0.0091 


13 


8 


0.1154 


0.0033 


0.0097 


16 


7 


0.1125 


0.0066 


0.0110 



electron beam [321 will make it possible to extend the range of Q 2 up to ~ 15 (GeV/c) 2 . In our 
opinion this will allow a drastic improvement in the quality of the experimental Nachtmann 
moments at Q 2 > 5 (GeV/c) 2 for n > 2. 

From the theoretical point of view the main improvements to be developed in our pro- 
cedure are: i) the consideration of more Nachtmann moments in our analysis in order to 
check the stability of the extracted values of a s (Mz) against the order of the moment; ii) 
the inclusion of all the NNLO effects with a consistent resummation of large n-logarithms 
at NNLL accuracy; and hi) the consideration of the full singlet and NS evolutions for the 
moments with n > 2. 

4 Conclusions 

An exploratory study for a new determination of the strong coupling constant a s {M z ) 
from existing world data on the structure function F2 of the proton in the Q 2 -range 5 < 
Q 2 (GeV/c) 2 < 120 has been presented. The main features of our approach are: 1) the use 
of low-order Nachtmann moments evaluated with a direct contribution from data larger than 
70% of the total; 2) the inclusion of high-order perturbative effects through the soft gluon 
resummation technique at next-to-leading-log accuracy; 3) a direct control over higher-twist 
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Table 5: The same as in Tabled but in case of the NLO approximation for the leading twist. 



moment 


a s (M z ) 


stat. 


syst. 


M 4 


0.1192 


0.0016 


0.0007 


M 6 


0.1221 


0.0013 


0.0020 


M 8 


0.1258 


0.0018 


0.0083 






average 


0.1221 


0.0009 


0.0031 


M 4 + M 6 


0.1209 


0.0010 


0.0015 



effects; and 4) the independence from any specific choice of the x-shape of the input parton 
distributions. 

At next-to-leading order we get a s {M z ) = 0.1209 ± 0.0010 {stat.) ±0.0015 (syst.) with a 
significant dependence upon the order of the moment used. Including soft gluon effects we 
obtain a s (M z ) = 0.1188 ± 0.0010 (stat.) ± 0.0014 (syst.) with a remarkable better stability 
against the order of the moment. Our findings compare positively with the total and DIS 
world averages, a s (M z ) = 0.1183 ± 0.0027 and a s (M z ) = 0.119 ± 0.002 (exp.) ± 0.003 (th.), 
respectively 

Thus, the new approach we are proposing appears to work quite well and it may represent 
an interesting tool for providing a precise determination of a s (M z ) from DIS data. Direc- 
tions for future improvements have been discussed. From the experimental side we need an 
improvement in the quality of the data, at least in the Q 2 range from rs 5 to w 20 (GeV/c) 2 , 
where the largest sensitivity to scaling violations is expected. The upgrade of Jefferson 
Lab to a 12 GeV electron beam and the use of the CLAS detector might provide such an 
improvement, offering the possibility to a very precise determination of the experimental 
Nachtmann moments of the proton (and deuteron) structure functions for n > 2. From the 
theoretical point of view the comparison of our findings with the recent NNLO results of 
Refs. |4J and indicates quite clearly the need of including all the NNLO effects with a 
consistent resummation of large n-logarithms at NNLL accuracy. Moreover, the full singlet 
and non-singlet evolutions should be considered for the moments with n > 2 as well as more 
Nachtmann moments should be included in our analysis in order to check the stability of 
the extracted values of a s (M z ) against the order n of the moment. Once such a stability is 
fully reached, the uncertainties on the extracted value of a s (M z ) can be further reduced by 
analyzing the scaling violations of all the moments with n > 2 simultaneously. 
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Figure 5: Values of a s (Mz) extracted in this work from the moments M n (Q 2 ) at n = 4, 6 and 
8. Full dots and open squares correspond to the results obtained using the SGR formula and 
within the NLO approximation for the leading twist, respectively. The vertical bars represent only 
the statistical uncertainties due to the fitting procedure. Note that the systematic errors are of the 
same order of magnitude as the statistical ones for n = 4 and n = 6, but much larger for n = 8 
(see Tables\^ and\$ . The solid line and the shaded area are the world average value and its global 
uncertainty, as given in Ref. J^. 
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